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' Abstract. We propose a non-equilibrium hadronic model to interpret the observed excitation func- 

O tion of kaon-to-pion ratios in nuclear collisions at AGS and SPS energies. The crucial assumption of 

our model is that due to stronger stopping at lower energies the lifetime of the fireball is prolonged 
' because the system has to build up the longitudinal expansion from internal pressure. 
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THE HORN 



> 

The excitation function of the ratio (K + ) / (n + ) which exhibits a sharp peak (a "horn") 
at projectile energies of 30 AGeV is certainly one of the most intriguing experimental 
results from the energy scan performed at the CERN SPS [1]. Statistical hadronisation 
model expects a maximum connected with the transition from baryon to meson domi- 
nated energy regime but cannot reproduce the observed sharpness of the peak |2D - The 
sharp maximum together with the excitation functions of (K~) / and (A) / (k) is 
not reproduced in hadronic transport codes either i3j,|4|]. The only successful, though 
schematic, interpretation of the data uses the framework of the so-called Statistical 
model of the early stage (SMES) [5]. The model is based on the assumption that the 
primordial particle production is realized according to grand-canonical equilibrium dis- 
tribution and the steep descent of the peak corresponds to transition to deconfined phase 
via a mixed phase. A kinetic description of the excitation function, also including phase 
transition in the region of the peak, has been proposed recently Jg]. 

Concerning the SMES, a question appears whether it is realistic to assume a chemical 
equilibrium thus early in the collision? Moreover, if the peak is to be regarded as 
a signature for deconfinement, are all hadronic scenarios safely excluded? We shall 
explore the possibility of reproducing the "horn" in a hadronic non-equilibrium model. 

Let us recall the two important quantities which regulate the final amount of produced 
strangeness. Firstly, it is the energy density, because the rates of reactions producing 
strangeness depend on the energy of particles in the incoming channel. Secondly, it is 
the total lifespan of the fireball since strange particles are out of chemical equilibrium 
and the relative strangeness content, most probably, grows with time. Thus we propose 
the following scenario potentially leading to the (K + ) / (n + ) "horn": at lower energies 
(< 30AGeV) the increase of the excitation function is due to an increase in the energy 
available for the strangeness production. With the further increase of the collision energy 



the nuclear stopping power decreases. Since we know from Mj -dependence of HBT 
radii that the longitudinal expansion pattern at freeze-out looks roughly the same at all 
energies, there must be accelerated longitudinal expansion and it must last longer at 
lower energies with stronger stopping. We put the sharp decrease of the ratio (on the 
right-hand side of the "horn") into connection with shorter lifetime of the fireball in 
more energetic collisions. The fireball just has less time to "cook up" strangeness and 
therefore the relative strangeness yield is lower. This hypothesis will be tested with a 
kinetic model. 



THE MODEL 

As we only want to calculate the ratios of total yields, and not the yields themselves, it 
is enough to study the evolution of the spatially averaged densities of individual species 
only. The spatial distribution is inessential here. We will in particular investigate the 
evolution of kaon density. Assuming that kaons are kept in thermal equilibrium with a 
fireball medium, the variation of the kaon density follows from the simple relation 

dn K _ d N K __ N K 1 dV 1 dNx 
~dx ~ dxV ~ ~^VV ~dx + V ~dx" 

where N% and V are the total number of kaons and the volume of the system, respectively, 
and T is the time in the co-moving frame. The second term on the right-hand side of this 
equation expresses the change of density due to chemical reactions. It can be split into 
gain term and loss term 
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In the gain term we sum over all processes producing a kaon, in the loss term these 
are the processes destroying a kaon. The relative-velocity-averaged cross-sections are 
multiplied with densities of the incoming species. 

The first term on the right-hand-sides of eqs. and © includes the expansion rate 
and stands for the density change due to expansion. In a simulation of the collision 
(hydrodynamic or transport), the expansion rate is obtained naturally. Here, we shall 
adopt a parametrisation of the expansion. Though in this way we do not make a direct 
contact to the underlying microscopic equation of state, it gives us a possibility to 
explore many different evolution scenarios and their impact on data. 

The ansatz for the energy density and the baryon density which will be used reads 

{Eoil-ax-bx 2 ) for x < z s 
£ o for x>x (3a) 

fps,o(l -ax-bx 2 ) 8 for x < x s 
PbW = U- for t>t <3b) 

In the above expressions the parameters can be tuned. Two of them are constrained by the 
requirements that the functions are continuous together with their first time derivatives. 



Initially (t < t s ) the fireball expands acceleratedly with a quadratic dependence on time. 
Then, at x > x s , the power-law expansion is dictated by the M^-dependence of HBT 
radii. Note the shift by To: it allows to delay this type of expansion. Such a delay would 
be unobservable in the Mj_ -dependence of i?i on g which is standardly used to argue for 
a short fireball lifetime 0]. The "lifetime" measured by ^i on g would correspond in our 
case to the difference (T tota i — To). Note also that the two parameterisations for energy 
and baryon density have similar forms and differ essentially by the exponent 8 which is 
given by the assumed equation of state. 

By tuning the parameters a and b we specify the initial rate of density decrease and 
the pace on which the decrease accelerates. This way we quantify the stopping and re- 
acceleration. We can construct a large class of models "between Landau and Bjorken 
scenarios". 

Densities of kaons (i.e. K + , K°, K* + , K*°) are calculated according to eq. ©. We 
include the following processes: 
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For those processes with two particles in final state also inverse reactions have been 
included which annihilate kaons. 

Species with negative strangeness must balance the total strangeness of the system to 
zero. Reactions which just swap the strange quark between them are quick. Therefore, 
we assume that these species are in chemical equilibrium with respect to each other, 
while the strangeness-weighted sum of their densities is given by the requirement of 
strangeness neutrality. 

All non-strange species are assumed to be chemically equilibrated. We assume no 
antibaryons at these energies. Their influence would be highest at the top SPS energy 
y/SNN = 17 AGeV where we expect 10% error due to their neglect. In a schematic model 
like ours this is acceptable. 

In a model set up in this way, the total produced amount of strangeness is mainly 
controlled by the lifetime (as we shall show) and slightly less by the temperature. The 
rate at which strange quarks are distributed among K and A is fixed by the final 
temperature. 

We try to choose the parameterisations © in such a way that the final state resulting 
from fits to chemical composition at different energies Jsh is reached. The initial amount 
of strangeness is estimated from pp, pn, and nn collisions H]. 



RESULTS 



We show an example of our results in Figure [T] Each point in that figure corresponds to 
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FIGURE 1. Calculated ratios {K + )/{lt + ), (K~) /{it'), (A) /{it) as functions of the total lifetime of the 
fireball. Different curves correspond to different initial energy densities. Calculation for Pb+Pb collisions 
at projectile energy of 30 AGeV. The bands show values accepted by data 1 1]. 



a different evolution scenario; they differ by total lifetime and initial energy density. 
Though the presented results are obtained for Pb+Pb collisions at projectile energy 
of 30 AGeV, it applies generally that dependence on lifetime is crucial, while the 
dependence on initial energy density is less important. The latter gains more weight 
at lower energies and is completely irrelevant at the highest SPS energy. It mainly stems 
from the temperature dependence of reactions including nucleons like % + N — > A + K 
which make up a bigger share of the total production rate at lower energies where the 
baryon density is higher. 

Summary of the values of total lifetime and initial energy density allowed by com- 
parison to data is presented in Figure El The more important limitations are put on the 
lifetime. In general, for SPS energies our hypothesis is confirmed that with increasing 
energy of the collision lifetime of the fireball decreases. Only at the lowest studied en- 
ergy, 1 1.6 AGeV (AGS), we miss the (A) / {n) ratio by about 3 standard deviations, see 
Figure |3j This is caused by much too high temperature in the final state, which can be 
fixed by re-parameterising the time dependence of the energy density. Otherwise, we 
reproduce the data well. 

CONCLUSIONS 

The excitation function of the ratios {K + ) / {n + ), {K~) / {n~) and (A) / {n) can be repro- 
duced in a non-equilibrium hadronic scenario. The crucial assumption is that the total 
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FIGURE 2. Values of initial energy density and total lifetime accepted by data QJ]. 
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FIGURE 3. Results obtained with parameter sets indicated in lower right panel of Figure [2] compared 
to data. 



lifetime of the fireball is a decreasing function of the collision energy. 

Thus we have made a specific prediction for the evolution dynamics of the fireball. 
Such a prediction must now be checked against other observables. Although we were 
lead by the knowledge of all particle abundances, transverse momentum spectra and 
HBT radii, a careful analysis of these data in the framework of the proposed model must 
be performed. Furthermore, dilepton spectra gained an unprecedented accuracy recently 
iTToL [Til H2I1 . They are sensitive to whole fireball history ifoll and could provide the 
crucial test of our hypothesis. 
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